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Heavy quarkonium production in hadron collisions pro- 
vides critical tests of quantum chromodynamics (QCD) 
because it involves both short-distance and long-distance 
contributions. The cross sections for direct charmonium 
and T production measured at the Tevatron in pp colli- 
sions [l| greatly exceeded the predictions from leading- 
order "color singlet" models [31. Further measurements 
of spin alignment in J/yj, "0' [3 an d Y production Q 
proved to be in dramatic disagreement with predictions 
of "color octet" models that were developed to explain 
the cross section results [5j. While spin alignment can 
provide very sensitive tests for QCD models, the mea- 
sured elements of the spin density matrix for the spin-1 
T states depend critically on the choice of coordinate 
frame. 

Recently, it has been pointed out that improved ex- 
perimental measurements are needed to clarify this situ- 
ation Q. In general, the angular distribution of the fi + 
in the rest frame of an T — > decay can be written 
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as 



7] 



dN 

In 



1 + Xg cos 2 9 + X v sin 2 9 cos 2ip + Xg v sin 29 cos if, 

(1) 

in which 9 is the polar angle measured with respect to 
a quantization axis, and <f is the azimuthal angle mea- 
sured with respect to the production plane containing the 
T and the beam axis. The coefficients are directly related 
to the elements of the spin-density matrix for the ensem- 
ble of T states observed 0. Previous studies at hadron 
colliders have measured only Xg in the s-channel helicity 
frame, where the quantization axis coincides with the di- 
rection of the T momentum, and cannot be transformed 
into different coordinate frames without additional in- 
formation. This precludes mo del- independent compar- 
isons of results obtained in different coordinate frames 
or experimental environments. Improving on this exper- 
imental situation requires not only measuring all three 
coefficients, but also carrying out these measurements in 
multiple coordinate frames. This would allow a com- 
parison of quantities such as A = (Xg + 3X v )/(l — X v ), 
which should have the same value in different coordi- 
nate frames 0] . Such a test would provide an important 
demonstration that observations have not been seriously 
biased due to poor determination of experimental accep- 
tance or subtraction of highly non-isotropic backgrounds. 

In this Letter, we report on the first analysis of angu- 
lar distributions of muons from T(1S,2S,3S) — > 
decays produced in pp collisions carried out using this 
formalism: the distributions are quantified in both the 
s-channel helicity frame and in the Collins-Soper frame, 
which approximates, on average, the direction of the ve- 
locity of the colliding beams. It is also the first analysis 
to provide information on the angular distributions of 
muons in decays of the T(35) state which is more likely 
to be produced directly, rather than as a decay product 
of higher mass quarkonium states. 

The T — > decays were collected using the CDF II 

detector, which reconstructs charged-particle tracks and 
measures their momenta using a six-layer silicon strip de- 
tector [Io| and a large- volume drift chamber 11 1, both 



with approximate cylindrical geometry and positioned in 
a 1.4 T solenoidal magnetic field. The tracking detectors 
are surrounded by calorimeters and three separate muon 
detector subsystems. The CMU system 12j consists of 
four layers of drift tubes that are located outside the 
hadron calorimeter and cover the central range of pseu- 
dorapidity \rj\ < 0.6. The CMP muon system 13j com- 
prises more layers of drift tubes and scintillator placed 
behind additional steel absorber material and covers ap- 
proximately \rj\ < 0.4. Extended muon coverage in the 
forward region 0.6 < \rj\ < 1 is provided by the CMX sub- 
detector which is also constructed from scintillator and 
drift tubes. 

A three-level online event selection system (trigger) is 
used to identify events that contain oppositely charged 



dimuon candidates. The level- 1 trigger requires two 
tracks with pt > 1.5 GeV/c to be identified in the track- 
ing chamber and to be geometrically correlated with ac- 
tivity in the CMU or CMX muon systems [3 . The level- 
2 trigger requires the muons to have opposite charge, and 
it requires that one of the muons have px > 3 GeV/c 
and that it be detected in both the CMU and CMP sys- 
tems. Events satisfying the level-2 trigger are passed to 
the level-3 trigger system, which employs a version of the 
full event reconstruction software optimized for speed. 
The event selection used in the level-3 trigger requires the 
presence of two oppositely charged muon candidates with 
invariant mass in the range 8 < m(/i + /i~) < 12 GeV/c 2 . 
It requires one to be reconstructed in both CMU and 
CMP systems with px > 4 GeV/c and the other to be 
reconstructed in either the CMU or CMX detectors with 
Pt > 3 GeV/c. In this Letter, we refer to the trigger 
scenario that selects two central muons as CC, and the 
scenario that selects one central and one forward muon as 
CF. From 6.7 fb _1 of integrated luminosity, the combina- 
tion of these triggers provides an event sample containing 
approximately 550 000 T(15), 150 000 T(2S), and 76 000 
T(35) decays. 

The criteria used to select dimuon candidates closely 
follow those previously used in Ref. 15]. Muon candi- 
dates are reconstructed from tracks in the drift cham- 
ber that extrapolate to a track segment reconstructed in 
at least one of the muon detector systems. Geometric 
restrictions are imposed on muon candidates to ensure 
that they are contained in regions of the detector with 
well-measured trigger and track reconstruction efficien- 
cies. Efficiencies for the level- 1 trigger and for these se- 
lection criteria are measured using the unbiased track in 
J/ip — > decays that were recorded using a single- 

muon trigger. This analysis also makes use of informa- 
tion from the CMP muon system in the level-2 trigger. 
The efficiency for selecting such muons is measured using 
samples of J/ip — > decays obtained using triggers 

that required information from only CMU or CMX. 

The angular distribution analysis is performed sepa- 
rately in each of the 12 ranges of dimuon mass shown in 
Fig. [U The angular distributions of T decays are ana- 
lyzed in eight ranges of Pr(T) from to 40 GeV/c and are 
restricted to the central region of rapidity |y(T)| < 0.6. 
For a given range of transverse momenta, the sample of 
dimuon candidates is divided into two subsamples ac- 
cording to whether one of the muons is reconstructed 
precisely using measurements from the silicon detector 
and its extrapolated trajectory misses the beam axis by 
a distance \do\ > 150 /im. Events with at least one muon 
satisfying this requirement are referred to as the "dis- 
placed" sample since they are consistent with the pres- 
ence of a long-lived parent particle, which is a charac- 
teristic feature of the dimuon background arising from 
semileptonic decays of heavy quarks. These criteria do 
not bias the angular distribution and, since the displaced 
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FIG. 1: Mass distribution of fi + fj,~ candidates in the prompt 
(solid line) and displaced (dashed line) samples, with the 
ranges of invariant mass used to select signal and sideband 
events indicated. Hatched regions indicate the three mass 
ranges containing the T(nS) signals. 

sample contains almost no T signal, it provides a good 
description of the dimuon background that remains in 
the complementary "prompt" sample. We verify that 
this is the case by comparing the angular distributions of 
prompt and displaced samples projected onto the cos 9 
and (p axes for mass ranges in the sidebands of the T 
signals. The level of agreement is quantified by comput- 
ing the Kolmogorov-Smirnov statistic for each distribu- 
tion, which suggests that any differences in their observed 
shapes are consistent with statistical fluctuations. 

The prompt sample contains most of the T — > 
signal, while a small fraction is retained in the displaced 
sample due to the 30 jirsx resolution of the do measure- 
ment (see Fig. [IJ . The fraction f p that is present in the 
prompt sample is measured using a simultaneous binned 
likelihood fit to the dimuon mass distributions of both 
prompt and displaced samples. The T signals are de- 
scribed by Gaussian functions with common widths, their 
mass splittings constrained to the known values [16| and 
their yields scaled by f p and 1 — f p in the prompt and 
displaced samples, respectively. The value of f p ranges 
from 96-99% depending on the T pt and whether the 
candidate was recorded with the CC or CF trigger sce- 
nario. 

A second fit is performed to measure the mass- 
dependent ratio r(m) of the prompt and displaced mass 
distributions. This is similar to the first fit, but uses 
only the sidebands in the prompt sample, < 
9 GeV/c 2 or m(/i+//-) > 10.5 GeV/c 2 , to avoid the 
need to model the T(nS) line shapes. The mass dis- 
tribution of the background in the displaced sample is 
parametrized by a gamma function at low-px and by an 
exponential function at higher pt- The mass distribu- 
tion of background in the prompt sample is accurately 
described by the mass distribution in the displaced sam- 




<t» (cleg) <|> (deg) 



FIG. 2: Comparisons of projected angular distributions mea- 
sured in the Collins-Soper frame for prompt (histograms) and 
displaced (error bars) dimuon samples in the low-mass (a,c) 
and high-mass (b,d) sidebands for events recorded with the 
CC trigger. Displaced samples are normalized to the number 
of events in the corresponding prompt samples. The p-value 
of the Kolmogorov-Smirnov statistic comparing each distri- 
bution is shown. 



pie multiplied by a scale factor that varies linearly with 
mass, r(m) = a + bm, with the coefficients a and b deter- 
mined from the fit. The value of the function r(mr(is))i 
evaluated at the T(1S) mass, varies between 1.8-3.9 over 
the range of dimuon pt considered. In the subsequent 
analysis of angular distributions, the value of this func- 
tion and its uncertainty, both evaluated at the center of 
each mass range containing the T signals, are used to 
impose a Gaussian constraint on the background yield in 
the prompt sample. 

The displaced sample provides a good description of 
the angular distribution of background in the prompt 
sample. This observation is consistent with all back- 
ground arising from semilcptonic decays of heavy quarks, 
or any small non-heavy flavor background component 
having the same angular distribution in prompt and dis- 
placed samples. We observe good agreement between the 
angular distributions in the prompt and displaced back- 
ground, outside the T signal mass regions, even though 
the angular distributions change rapidly with dimuon 
mass and px- A typical example illustrating this com- 
parison is shown in Fig. [5] We then proceed to use the 
displaced sample to constrain the angular distribution of 
the background when analyzing the angular distribution 
of muons from T decays. 

In mass ranges containing the T signals, we perform a 
third simultaneous fit, in both the prompt and displaced 
samples, to the distributions of angles (cos 9, ip) collected 
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in 20 x 36 discrete intervals. Separate fits are performed 
for the two cases in which the angles represent the direc- 
tion of the positive muon with respect to the axes of the 
s-channel helicity frame or those of the Collins-Soper 
frame. For both frames, the angular distributions are 
described by linear combinations of probability density 
functions for signal and background components. We fac- 
tor these functions into an underlying angular distribu- 
tion, with that of the T signal parametrized using Eq. (flj , 
and an acceptance function that accounts for the geome- 
try of the muon detectors and the kinematic restrictions 
imposed by the trigger. The parameters in the under- 
lying angular distributions for the T signal and for the 
background are determined using a simultaneous, binned 
likelihood fit in which the expected numbers of events in 
each discrete angular interval are expressed as 

NdS p Ab(cos9i,<pj) ■ iu&(cos#j, tpj; A 6 ), (2) 
^± rv, jVr(l - f p )Ar (cos 9i,ifj) ■ w r (cos di,ipj] A T ) + 

N d Ab(cos6i,ipj) ■ w b (cos6i,cpj;\ b ). (3) 




cose <|> (deg) 



FIG. 3: Examples of projections of angular variables mea- 
sured in the Collins-Soper (a,b) and the s-channel helicity 
(c,d) frames for the range of invariant mass containing the 
T(15) signal 4 < p T < 6 GeV/c recorded with the CC trig- 
ger. The prompt sample is shown in black while the displaced 
sample is shown in grey. The distributions from collision data 
are shown with error bars and the results of the fits are indi- 
cated by the histograms. 



In these expressions, Nr and Nd are the T and displaced 
background event yields, f p is the fraction of the T sig- 
nal in the prompt sample, and s p is the ratio of the 
background yields in the prompt and displaced samples, 
which is Gaussian-constrained to r(m). The acceptance 
for signal At and background Ab, which are calculated 
separately, are described below. The modeling of the 
background angular distribution is improved by imposing 
the additional kinematic restriction \pt(^ + ) ~~ Pt(^)\ < 
(Pt(a* + aO — 0.5 GeV/c), which removes back-to-back 
muons that have large values of cos 9 in the s-channel 
helicity frame. This has a negligible effect on the T(nS) 
acceptance for pr(T) > 6 GeV/c. The angular distribu- 
tions for T signal wr and background Wb are described 
by sets of parameters Ax and A&. For the signal, wx has 
the form of Eq. (TTJ), whereas u>b(cos 9, <p; A&) has, in addi- 
tion to the terms in Eq. ([1]), an empirical A4 cos 4 9 term 
that parametrizes the background shape more accurately 
in some ranges of pt and invariant mass. The free pa- 
rameters in the fit are Ny, N^, s p , Ax, and A&, while f p 
is fixed to the value determined previously from the first 
fit to the dimuon mass distributions. 

The detector acceptance is calculated using a Monte 
Carlo (MC) simulation in which dimuon events gener- 
ated with isotropic distributions of decay angles are pro- 
cessed using the standard CDF II detector simulation 
and event reconstruction programs. Separate samples 
are simulated at fixed masses to calculate the acceptance 
for the three T(nS) signals, while the acceptance for the 
dimuon background is calculated using a continuum of 
dimuon invariant masses ranging from 8 to 12 GeV/c 2 . 



Figure [3] shows distributions for the data and the best 
fit model in one of the ranges of dimuon px in the mass 
range containing the T(1S) signal. The shapes of the 
projected distributions are primarily determined by the 
acceptance since the resolution with which cos 9 and (p 
are measured is much smaller than the binning chosen 
for the fit. However, the agreement between the data 
and the model in regions where the acceptance changes 
rapidly, demonstrates that resolution effects are modeled 
accurately in the detector simulation. No significant dis- 
crepancy between the data and the fit model is observed 
in x 2 tests applied to one dimensional projected distri- 
butions over the analyzed range of dimuon mass and px ■ 

Systematic uncertainties on the parameters Xg , X v and 
Xg v due to the limited precision with which the trig- 
ger and reconstruction efficiencies are determined are 
evaluated by repeating the analysis with acceptances re- 
calculated with all efficiencies simultaneously varied by 
±1(7. The resulting change in the fitted parameters pro- 
vides a conservative estimate of the sensitivity to the 
measured acceptance. Because the measured parame- 
ters depend on the estimated background in the prompt 
sample, an alternate, quadratic parametrization of the 
function r(m) was also investigated. The resulting small 
variations in the fitted angular distribution parameters 
are treated as an additional systematic uncertainty. Fi- 
nally, the contribution of the uncertainty in the fitted 
parameters due to the finite MC sample size used to cal- 
culate the acceptance was estimated using ensembles of 
MC simulations with the same size used in the analysis of 
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FIG. 4: Rotational invariant A as a function of pt(T) for 
the T(IS'), T(2S), and T(35') states. Values of A calculated 
in the Collins-Soper frame are indicated by dark lines, while 
those calculated in the s-channel helicity frame are indicated 
by grey lines and are horizontally offset to slightly larger pr 
values for clarity. 



the data. The uncertainties due to finite MC sample size 
and the determination of efficiencies are at most 30% of 
the size of the statistical uncertainty for the three T(nS) 
states, while the uncertainty due to the treatment of the 
prompt scale factor function is no more than 20% of the 
statistical uncertainty. 

Figure [4] shows the rotational invariant A which is cal- 
culated from the measured values of \g and X^ in each 
Pt range for both the Collins-Soper and s-channel he- 
licity frames. Uncertainties in A measured in the two 
coordinate frames are highly correlated. Monte Carlo 
simulations are used to calculate the expected sizes of 
differences between the two values of A and in most cases, 
the observed deviations are found to be consistent with 
purely statistical fluctuations. A systematic uncertainty 
derived from the difference between A measured in the 
two coordinate frames is only significant for the lowest 
three pr ranges of the T(3S*). In the lowest pt range, 
the values of A measured for the T(3S) differ by 2.4 er, 
without accounting for systematic uncertainties, and is 
the only case where the angular distribution is observed 
to be significantly non-isotropic. However, we can find 
no evidence to suggest that this is due to a bias or sys- 
tematic effect since we do not see a similar trend for the 
IS and 25 states and can find no anomalous behavior in 
any of the underlying distributions. 

The values of A w suggest that the decays of all 
three T(nS) resonances are consistent with an unpolar- 
ized mixture of states. Table Q] lists the values of Ag 
measured in the s-channel helicity frame for the T(15), 
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FIG. 5: Comparison of the Xe parameter measured for T(15*) 
decays in the s-channel helicity frame (solid symbols) with 
previous measurements from the CDF [J] (open circles) and 
the DO [18I ] (open triangles) experiments. 



T(25), and T(35) states, with the systematic uncertain- 
ties described above added in quadrature [l7|. Figure [5] 
shows a comparison of the Xg parameter, measured for 
the T(IS') state in the s-channel helicity frame, with 
previous measurements. The current result is found to 
be statistically consistent with the previous measure- 
ment from CDF which was made for \y\ < 0.4 at 
y/s = 1.8 TeV rather than \y\ < 0.6 and Vs = 1.96 TeV. 
Restricting the current measurement to \y\ < 0.4 does not 
change the results appreciably. The current T(15) result 
is inconsistent with the previous measurement from the 
DO experiment [3| at the level of 4.5a. 

In conclusion, we have measured the angular distri- 
butions of muons from T(l-S), T(2S), and T(35) decays 
with \y\ < 0.6 and in several ranges of transverse momen- 
tum up to 40 GeV/c. We find that the decay-angle distri- 
butions of all three T(nS) states are nearly isotropic, as 
was suggested by previous measurements [4( in the case 
of the T(1S). This is the first measurement to simulta- 
neously determine the three parameters needed to fully 
quantify the angular distribution of T(nS) — > fJ. + fJ>~ de- 
cays. This is also the first analysis to present information 
on the angular distribution of T(3S') mesons produced in 
high energy pp collisions. 
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pr (GeV/c) 


Xg (T(15)) 


Xg (T(25)) 


Xg (T(35)) 


0-2 


0.01 ± 0.10 ± 0.05 


-0.07 ±0.20 ±0.05 


0.32 ± 0.36 ± 0.12 


2-4 


-0.08 ±0.04 ±0.03 


-0.29 ± 0.10 ± 0.06 


-0.19 ±0.17 ±0.08 


4-6 


-0.11 ± 0.04 ± 0.03 


-0.11 ± 0.10 ± 0.05 


-0.01 ± 0.19 ± 0.09 


6-8 


-0.12 ±0.06 ±0.04 


-0.09 ±0.12 ±0.07 


-0.48 ±0.14 ±0.13 


8-12 


-0.15 ±0.06 ±0.02 


-0.18 ±0.11 ±0.03 


0.07 ±0.18 ±0.04 


12-17 


-0.13 ±0.08 ±0.04 


-0.09 ±0.15 ±0.08 


0.14 ±0.22 ±0.07 


17-23 


-0.23 ±0.13 ±0.08 


0.06 ±0.25 ±0.12 


0-i4±g;H ±0.12 


23-40 


-0.21 ±0.24 ±0.11 


O.OOtg-j^ ±0.16 


0.14l° j° ± 0.13 



TABLE I: Values of Xg, measured in the s-channel helicity frame for each range of pt(T). The first uncertainty is statistical 
while the second uncertainty is systematic and is described in the text. 
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